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Abstract 

When a suspension consisting of dielectric particles having nonlinear charac- 
teristics is subjected to a sinusoidal (ac) field, the electrical response will in 
general consist of ac fields at frequencies of the higher-order harmonics. These 
ac responses will also be anisotropic. In this work, a self-consistent formalism 
has been employed to compute the induced dipole moment for suspensions in 
which the suspended particles have nonlinear characteristics, in an attempt 
to investigate the anisotropy in the ac response. The results showed that the 
harmonics of the induced dipole moment and the local electric field are both 
increased as the anisotropy increases for the longitudinal field case, while the 
harmonics are decreased as the anisotropy increases for the transverse field 
case. These results are qualitatively understood with the spectral represen- 
tation. Thus, by measuring the ac responses both parallel and perpendicular 
to the uniaxial anisotropic axis of the field- induced structures, it is possible 
to perform a real-time monitoring of the field-induced aggregation process. 

PACS Number(s): 83.80.Gv, 72.20.Ht, 82.70.Dd, 41.20.-q 
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I. INTRODUCTION 



When a suspension of highly polarizable dielectric particles is subjected to an intense 
electric field, the induced dipole moments cause the particles to form chains along the field 
direction, resulting in complex anisotropic structures. These chain structures alter the ap- 
parent viscosity of the suspension by several orders of magnitude, and are the basis of the 
electrorheological (ER) effect. Under the action of the ER effect, the arrangement of parti- 
cles, or the shape of the lattice changes and the local electric field deviates from the Lorentz 
cavity field [|T|. As a result, these structures have anisotropic physical properties, such as in 
the effective conductivity and permittivity [§], and in the optical nonlinerity enhancement 
1^. The rapid field-induced aggregation and the large anistropy in their properties render 
this field-induced-structured material potentially important for technological applications 

i- 

Moreover, the applied electric field used in most ER experiments is quite high, and 
important data on nonlinear ER effects induced by a strong electric field have been reported 
recently 0. A convenient method of probing the nonlinear characteristics of the field- 
induced structures is to measure the harmonics of the nonlinear polarization under the 
application of a sinusoidal (ac) electric field. Thus by measuring the ac responses both 
parallel and perpendicular to the uniaxial anisotropic axis of the structures, it is possible to 
perform a real-time monitoring of the field-induced aggregation process. 

Recently, the effect of nonlinear characteristics on the interparticle force has been an- 
alyzed in an ER suspension of nonlinear particles [§] and further extended to a nonlinear 
host medium 0. The nonlinear characteristics are due to the field-dependence of the di- 
electric constants of the materials used in ER fiuids, which have a constitutive relation 
D = eE + xE^. When a nonlinear composite with nonlinear dielectric particles embedded 
in a host medium, or with a nonlinear host medium is subjected to a sinusoidal field, the 
electrical response in the composite will in general be a superposition of many sinusoidal 
functions . It is natural to investigate the effects of a nonlinear characteristics on the ac 
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response in an ER fluid which can be regarded as a nonhnear composite medium . The 
strength of the nonhnear polarization should be reflected in the magnitude of the harmonics. 



In this work, we will develop a perturbation expansion method and a self-consistent 



theory |]T2[ to calculate the ac response of a nonlinear compsite. The paper is organized 
as follows. In the next section, we calculate the dipole moment of the polarized spherical 
particles in an anisotropic composite and extract the harmonic response. In section III, 
we perform a series expansion of the local fleld inside the particles from the self-consistent 
solution, in an attempt to obtain analytic expressions of the higher harmonics of the dipole 
moment. Numerical results are performed in section IV to validate the analytic results. 
Discussions on the results will be given. 

II. NONLINEAR POLARIZATION AND ITS HIGHER HARMONICS 

We flrst examine the effect of a nonlinear characteristics on the induced dipole mo- 
ment. We concentrate on the case where the suspended particles have a nonlinear dielectric 
constant, while the host medium has a linear dielectric constant 62. The nonlinear charac- 
teristics gives rise to a fleld-dependent dielectric coefficient ]13|. In which case, the electric 



displacement-electric fleld relation inside the spheres is given by: 

Di = eiEi + x(^i)Ei = eiEi, (1) 

where ei and x ^"^^ the linear coefficient and the nonlinear coefficient of the suspended 
particles respectively. 

This constitutes an approximation: the local fleld inside the particles is assumed to be 
uniform and the assumption is called the decoupling approximation ||12[ . It has been shown 
that such an approximation yields a lower bound for the accurate result for the local fleld 
|T^. We further assumed that both e and x ^'^^ independent of frequency, which is a valid 
assumption for low- frequency processes in ER fluids. As a result, the induced dipole moment 
under an applied fleld 'Ei = E{t)z. is given by: 
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P = eea'bE{t), (2) 



where b is the field-dependent dipolar factor and is given by: 

~ ei - £2 ei + x{Ef) - €2 



(3) 



ei + 2e2 + x{E!) + 262' 
In order to obtain the effective dielectric constant ie, we invoke the Maxwell-Garnett 
approximation (MGA) for anisotropic composites For the longitudinal field case when 
the ac field is applied along the uniaxial anisotropic axis, the MGA takes on the form 

^--^^ = flLZ^ U) 

PL~ee + i^-PL)e2 ei + 2e2' 

whereas for an transverse field case when the ac field is applied perpendicular to the uniaxial 

anisotropic axis, the MGA reads 

~ €2 _ J Cl ~ ^2 

/?Tee + (3-/3r)e2 ~^ei + 2e2' ^ 
where jSi and /St denote the local field factors parallel and perpendicular to the uniaxial 
anisotropic axis respectively, and / is the volume fraction of particles. These factors are 
defined as the ratio of the local field in the particles to the Lorentz cavity field |l|]. For 
isotropic composites, /?l = /3t = 1, while both and Pt will deviate from unity for an 
anisotropic distribution of particles in composites. The f3 factors have been evaluated in 
a tetragonal lattice of dipole moments |I| and similar quantities were calculated in various 
field- induced-structured composites |T^; they satisfy the sum rule: 

/3l + 2/3t = 3. 

In what follows, we denote /Sl as f3 for notation convenience. 

On the other hand, the effective nonlinear dielectric constant can be given by [|T^ 

/<~1 , (l-/)f2 



E^{t)V Jv ' " ' ' " E^it)' " E^(t) 

where V the volume of the composite and / is the volume fraction. Accordingly, the local 
electric field inside the spheres can be expressed in terms of the derivative of le with respect 
to ii'. 



(^1^) = (7) 

If we apply a sinusoidal electric field, i.e. E(t) = EQsm{ut), the induced dipole mo- 
ment p will depend on time sinusoidally, too. By virtue of the inversion symmetry, p is a 
superposition of odd-order harmonics such that 

p = Pcj sin ut + sin Sut H . (8) 

Also, the local electric field contains similar harmonics 

^/{Efj = sin ujt + E3^sm3ujt + ---. (9) 

These harmonic coefficients can be extracted from the time dependence of the solution of p 
and Ei{t). 



III. ANALYTIC SOLUTIONS 



In what follows, we will apply two methods to extract the harmonics of the induced dipole 
moment and the eletric field: the perturbation expansion method (PEM) [|ri|,|l^ and the 



self-consistent (SC) theory ||T2[. The self-consistent theory can deal with the case of strong 
nonlinearity, while the perturbation expansion method is applicable to weak nonlinearity 
only, i.e. x(-^i) I5 limited by the convergence of the series expansion. In the case of 
weakly nonlinear response, the SC theory should agree with PEM. 

A. Perturbation expansion method 



We expand p and JxiEf) into a Taylor expansion 



P = a'E{t)Y.a^{x{Ef)r, (10) 

n=0 



m=0 



X{EI) = ^^j^ E d^{x{El)r, (11) 



5 



where 



1 9" ^ ~ 1 d"^ Idee 



For weak nonlinearity, we can rewrite Eqs.(lO) and (11), keeping the lowest orders of 



- ^^^(do + d,x{El)l (12) 
P = a^aoE{t) + ja^aidlxE^t) = hiE{t) + h^xE^t). (13) 

In the case of a sinusoidal field, we can expand E^{t) in terms of the first and the third 

harmonics. The comparison with Eq.(8) yields the harmonics of the induced dipole moment: 

3 1 

^hiEo + -hsxEo, P3uj = --hsxEo- 

Similarly, we find the harmonics of the local electric field 

X'^'E^ = nx"'E, + \j,{x"'E,)\ x"'E,^ = -i,3(x^/^i?o)^ 

where 

In the above analysis, we have used the identity sin^ ujt — ^ sin ujt — \ sin 2>ujt to obtain 
the first and the third harmonics. Similar analysis can be used to extract the higher-order 
harmonics, by retaining more terms in the series expansion. 

B. Self-consistent theory 

Eq.(7) is actually a self-consistent equation since its right-hand side depends on the local 

field itself. It must be solved sclf-consistently, and the local electric field reads: 

r-™- ^ -2{2,fl^R + 2V3(9gv/T + ^/?,^AK-^ + 27Q^T)V^ 

1/ 62/3VT(9g\/T + \/3x/4R^T27Q^)V3 ^ ' 

where 

g = 3e2x'/'£^(i), i? = eiT + e2(2 + //3), T = 1 - 
For a sinusoidal applied electric field, we determine numerically the harmonics of the 
induced dipole moment and the local electric field, which is a subject for the next section. 
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IV. NUMERICAL RESULTS 



In this section, we perform numerical calculations to investigate the effects of a nonlinear 
characteristics on the harmonics of the induced dipole moment and the local electric field. 
As shown in Section III, the harmonics of the induced dipole moment and the local electric 
field are both related to the field-induced anisotropy parameter (3. In order to investigate the 
effect of field-induced anisotropy on the harmonics, we perform the numerical calculations. 
Without loss of generality, we let / = 0.09, a = 1, ei = 10, and e2 = 1 for model calculations. 
We let po = ^eO'^bEQ to normalize the induced dipole moment. 

In Fig.l, we plot the normalized harmonics versus (3 by using the SC theory for the 
longitudinal electric field. It is evident that, for increasing /3, the harmonics of the induced 
dipole moment and the local electric field are both increased. Similarly, in Fig.2 the same 
quantities are plotted but for the transverse electric field. We find that the harmonics are 
decreased as P increases. 

Both Fig.l and Fig.2 also show that the harmonics are strongly dependent on the nonlin- 
ear response of the suspended particles. Moreover, increasing the nonlinear response leads 



to an increase in the harmonics. This is consistent with the results of our recent work [jT5 
in which a pair of nonlinear particles suspended in a linear host was considered. 

In Figs. 3 and 4, we compare the SC theory and PEM. For weak nonlinearity, xEq = 0.09, 
the PEM results coincide with the SC results, while for moderate nonlinearity x^q = 0.9, 
there are some deviations. In both cases, we find that the PEM results are in good agreement 
with the SC results for weak nonlinearity. For strong nonlinearity, there are large deviations 
(not shown here) and the SC theory must be used. 

The f3 dependence of the local field can qualitatively be understood with the spectral 
representation [jT6[. From Eq.(7), denoting s = (1 — ei/e2)^^, we find that Ei = sEo/{s — si), 
where si = (1 — (3f)/3. Using the numerical values of our model calculations, we find 
s = —1/9, and Ei = (l/9)-E'o/(l/9 + si). When /? increases, si decreases, leading to an 
increase in the local electric field Ei. 
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DISCUSSION AND CONCLUSION 



Here a few comments on our results are in order. In the present study, we have examined 
the case of nonhnear particles suspending in a linear host. We may extend our considerations 
to a nonlinear host medium ||^. In this case, preliminary results show that qualitatively 
similar yet more complex behaviors in the ac response have been observed. 

So far, we have not considered the frequency dependence of the particle dielectric con- 
stant. In a realistic situation, the dielectric constant of the particles can decrease with the 
increase of the frequency. For simplicity, we may adopt the Debye relaxation expression for 
ei. Preliminary results show that the ratio of the third to first harmonic decreases with 
frequency, results that are in accord with recent experimental data of Ref. 

To our knowledge, an accurate evaluation of the local field factor (3 during field-induced 
aggregation is lacking. We suggest that the Ewald method can be extended to compute 
the (3 factor in computer simulation of ER fluids. This is a formidable task for the future. Our 
present theory may be applied to the electrorheology of milk chocolate [|1^. By measuring 
the nonlinear ac response, it may be possible to monitor the food production processes. 

In conclusion, a self-consistent formalism has been employed to compute the induced 
dipole moment for suspensions in which the suspended particles have a nonlinear character- 
istics, in an attempt to investigate the anisotropy in the ac response. 
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FIGURES 

FIG. 1. The harmonic responses of the induced dipole moment and the local electric field versus 

P for the longitudinal field case: XE^ = 9 (dotted fines), x^o = 25 (solid lines) and x^o = 36 
(dashed lines). 

FIG. 2. Same as Fig.l, but for the transverse field case. 

FIG. 3. The comparison between the SC theory and PEM with xEq = 0.09 and x-Eq = 0.9 for 
the longitudinal field case. 

FIG. 4. Same as Fig.3, but for the transverse field case. 
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